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We study a surface model with a self-avoiding (SA) interaction using the canonical 
Monte Carlo simulation technique on fixed-connectivity (FC) triangulated lattices of 
sphere topology. The model is defined by an area energy, a deficit angle energy, and the 
SA potential. A pressure term is also included in the Hamiltonian. The volume enclosed 
by the surface is well defined because of the self-avoidance. We focus on whether or 
not the interaction influences the phase structure of the FC model under two different 
conditions of pressure Ap; zero and small negative. The results are compared with the 
previous results of the self-intersecting model, which has a rich variety of phases; the 
smooth spherical phase, the tubular phase, the linear phase, and the collapsed phase. 
We find that the influence of the SA interaction on the multitude of phases is almost 
negligible except for the evidence that no crumpled surface appears under Ap = at least 
even in the limit of zero bending rigidity a — ¥ 0. The Hausdorff dimension is obtained in 
the limit of ct— >0 and compared with previous results of SA models, which are different 
from the one in this paper. 
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1. Introduction 

Over the past few decades, a considerable number of studies have been conducted 
on the surf ace models. The model was constructed for strings and membranes 
| 1 | 2 | 3 | 4 | 5 |617| anc [ jj. wag jeflngj on t ne Das i s of the differential geometric notion of 
curvatures [Sllinrn. The so-called crumpling transition is a shape transformation 
between the smooth phase at sufficiently large bending rigidity a and the c ollapsed 
phase at a— >0, and it has l ong been studied both theoretically I 12 | 13 | 14 LL5| ancl nu- 
merically I 16 | 17 | 18 | 19 | 20 | 21 ] wh ne t ne transition is considered as a continuous one 
| 14 | 15 |^ a possibility that it is of first-order is p ointed out and renormalization 
group studies ^2 anc j recent numerical studies I 23 | 24 | 25 |26J pj-g^igt that the transi- 
tion is of first order. The tra nsiti on was observed in the canonical surface model on 
relatively large sized surfaces ^3. 

In addition to the smooth and the collapsed phases, a variety of phases including 
the tubular phase are observed in surface models ^3, which are defined by a one- 
dimensional bending energy on the cytoskeletal structure. A planar phase and an 
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oblong linear phase can be seen in a model , which is defined by a one-dimensional 
bending energy and the Nambu-Goto area energy on the fixed-connectivity (FC) 
surface. A surface model defined by a deficit angle energy also has a rich variety of 
phases including a tubular phase It must be noted that these phase transitions 
can be observed on relatively smaller surfaces in contrast to the above mentioned 
crumpling transition of the canonical surface model. 

To construct a surface model, the self-avoiding (SA) property should be taken 
into account if we focus on membranes ^ l 32 | 33 | 34 | 35 | 36 | 37 | 38 | 39 | 40 | 41|42 | However, 

numerical studies of the SA surfaces are very time consuming because of the non- 
local pr operty of the interactions. The simulations on such large surfaces like those 
in Refs. I 25 | 26 | are no t feasible on currently available computers. Nevertheless, 
the numerical studies on the above mentioned variety of phases in those specific 
models are considered to be feasible on the SA surfaces. 

Therefore, it is very interesting to study whether or not the SA interaction 
influences the phase structure in those models without SA interactions (phantom 
surface models). It is possible that the multitude of phases is strongly influenced 
by the SA i nteractions . In fact , no completely-collapsed phase is observed in FC 
SA surfaces | 33|34 | 35 | 36 | 37|38 | 42 | Moreover, the SA interaction is expected to play 
a non-trivial role in the membrane morphology even at the smooth phase. It was 
recently reported that SA property is essential for a variety of shapes of the so-called 
excess cone at high bending regime 

In this paper, we study the surface model in Ref. with a SA interaction 
on FC triangulated surfaces by using the canonical Monte Carlo (MC) simulation 
technique. The smooth spherical phase, the tubular phase, the line ar phase, and 
the collapsed phase are seen in the FC phantom surface model Our interests 
are focused on whether or not such a variety of phases, including the collapsed 
phase, are influenced by the SA interaction. Two different values of pressure Ap are 
assumed such that Ap is zero and small negative. 

This paper is organized as follows: in Section [2j we make a brief outline of 
the current results of the numerical studies of phantom surface models and SA 
surface models on triangulated surfaces. In Section|3l we define the model with a SA 
interaction, which is slightly different from the currently well-known SA interactions 
for numerical studies. The Monte Carlo simulation technique is shown in Section HI 
and the numerical results are presented in Section [5l We summarize the results in 
the final Section El 

2. Triangulated surface models 

2.1. Phase structure of phantom surface models 

In this subsection, we give a brief outline of the phantom surface models on tri- 
angulated lattices in R 3 and the current numerical results. We start with the con- 
tinuous model, which is given by the continuous Hamiltonian S = S± + aS% , where 
Si = J ^/gd 2 xg ab d a X^d b X^ and S 2 = (1/2) / y /gd 2 xg ab d a n fl d b n f " . Si is just iden- 
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tical with the action of Polyakov string l ^ l iU l, where X^ denotes a mapping from 
a two-dimensional surface M to R 3 and represents the surface position in R 3 , g ab 
is the inverse of the metric tensor g a \, of M, and g is the determinant of g a b- The 
variables (ari, X2) represent a local coordinate of M. The image X(M)(c R 3 ) is the 
surface, which is triangulated in numerical studies. The symbol in S2 denotes 
a unit normal vector of X(M), and £2 is called the bending energy, and a is the 
bending rigidity. 

If g a b is fixed to the Euclidean metric S a b and X(M) is triangulated by piece- 
wise linear triangles, then we have S — Si+aS 2 , Si = ^ (Ij-lj) , S2 — 
Y^ij (1 — Hi • rij), where Xi(£ R 3 ) in Si is the position of the vertex i, in S2 
is a unit normal vector of the triangle i. The FC model is statistical mechanically 
defined by the partition function 

Z Rx = I HdXiexp[-S(X)], (fixed), (1) 

i—l 

where the prime in J' YiiLi dXi denotes that the three-dimensional multiple inte- 
grations are performed by fixing the center of mass of the surface at the origin of R 3 
to remove the translational zero mode. We call the FC model defined by the ener- 
gies Si and S2 as the "canonical" surface model. It was reported that the canonical 
model on surfaces of sphere topology undergoes a first-order transition at finite a c 
between the smooth phase at a — ¥ 00 and the collapsed phase at a — > The 

role of the Gaussian bond potential Si is to make the mean bond leng th constant 

1901 

and, hence, Si can b e re placed by a Lennard- Jones type potential and also by 
a hard-wall potential ^Sl. In a surface model on triangulated lattices of the seminal 
paper Ref. US f Kantor and Nelson, Si is given by a hard-core and hard- wall po- 
tential. This type of potential can be used as a SA potential, which is described in 
the following sebsection. 

A variation of the canonical model is obtained by replacing Si with the Nambu- 
Goto area energy Sa — J2a ^-a, where A a is the area of the triangle A. Sa is also 
obtained from the above mentioned continuous Hamiltonian Si by fixing g a b as the 
induced metric gab — daX^dbX^ of the mapping X. We call a model as the Nambu- 
Goto surface model if the Hamiltonian includes Sa as the bond potential term. 
It is well-known that the Nambu-Goto model with the canonical bending energy 
S2 = J2ij ( 1 — n i • n i ) is ill-defined in the sense that no equilibrium configuration is 
obtained in the numerical simulations SH. The ill-definedness comes from the fact 
that the area A a is totally independent of the shape of A, and the oblong and 
very thin triangles, which are considered as singular triangles, dominate the surface 
configurations in the whole range of a. However, if the canonical bending energy S2 
is replace d by a deficit angle energy such as S 2 nt — J2i (^» — 00 ) or S 2 nt = l^i~0o| 
| 45 | 46 | 47 |^ ^ ne moc [ e i turns to be well defined except in the limit of a — > q 129] 
symbol Si in S 2 nt is the sum of internal angles of triangles meeting at the vertex i, 
and (f>o is a constant and fixed to (f>o = 27r if the surface is closed. The deficit angle 
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energy SJf' = — J2i l°S (^»/27r) is possible on closed surfaces such as a sphere I2S|48| 
Those deficit angle energies are called as the intrinsic curva ture energy. The reason 
of the variety of phases in the Nambu-Goto model with S™*^^ seems that both Sa 
and S™* are insensitive to the surface shape. In fact, Sa is independent of whether 
or not the surface is composed of almost-regular triangles or oblong triangles. S™* is 
also independent of whether the surface is planar or cylindrical. Nevertheless, both 
of the smoothphase and the collapsed phase are stable on the disk surface and 
on the torus ^4 In this, paper we study the Nambu-Goto surface model with the 
intrinsic curvature S^. 

We should comment on the reason why we use S™* — — log (<5i/27r) as the 
intrinsic curvature energy. The origin of S™* = — J2i 1°S (^»/2tt) is the measure 
factor qf in the integrations J YiiLi dXiqf in Z, where qt is the coordination 
number of the vertex i and <r(= 3/2) is a constant ED. By identifying qi with Si 
and extending the constant a to the variable coefficient a, we have the expression 
—a log (Si). Including the normalization factor 2tt, we have the curvature energy 

4 nt = -E«i°g(<V2^)- 

We should also comment on the fact that the mean value of Si is constant 
such that (Si /N) =3/2 even in the limit of b -> 0. The reason for (Si/N) = 
3/2 is understood from the scale invariant property of the partition function^. 
In fact, by rescaling the integration variable in Z such that X — > XX, we obtain 
Z{\) = A 3 ^- 1 ) f Yl? =1 dXiexp[-S(\X)}, where S(XX) = X 2 S 1 +bS 2 - The scale 
invariance of Z indicates that Z(X) is independent of A and, therefore, is represented 
by dZ{\)/d\\ x= i=Q. Thus, we have (Si/N)=3(N-l)/2N~3/2. 

A variety of phases can also be seen in a model, which is obtained by replacing 
Si and S2 of the canonical model with Sa and the one-dimensional bending energy 
S^, respectively^. In this case, S^ is sensitive to the surface shape, while Sa 
is not as mentioned above. 

A variation of the canonical model i s obtained also b y including fluidity, which 
represents a lateral diffusion of vertices I 6 | 16 | 17 | 19 | 20 | 21 | two-dimensional flu- 
idity is defined on dynamically triangulated surfaces, where the triangulation T 
is considered as a dynamical variable of the model. The partition function of the 
model with fluidity is thus given by 

,1 n 

Ztin =Y, Y[dX t exp [-S(X, T)] , (fluid) (2) 

r t=i 

where S(X, T) represents that S is dependent on the variables X and T, and 'Y^-j 
represents the sum over all possible triangulations. In the fluid model corresponding 
to the canonical model, we cannot see the transition, which is seen in the canonical 
model on FC spherical surfaces. This is expected from the phase structure of com- 
partmentalized fluid surfaces ^21, where the lateral diffusion is allowed only inside 
the compartment, which is a sublattice structure on the surface. In this compart- 
mentalized model, a first-order transition, which is considered to be identical to the 
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one in the canonical FC model, disappears if the compartment size Lc is increased. 
The homogeneous fluid surface is obtained from the compartmentalized fluid sur- 
face by maximizing Lc such that the surface is composed of a single compartment 
or the surface has no compartment. Thus, we understand that the transition, which 
is observed on the compartmentalized fluid surfaces at relatively small Lc, cannot 
be observed on the homogeneous fluid surfaces. The Nambu-Goto model with the 
intrinsic curvature energy is well-defined even on the fluid surfaces and has a variety 

loni 

of phases 1 ^. 

By combining two different sets of ball-spring systems, Boal and Seifert intro- 
duced a fluid surface model with cytoskeletal structures, which is a two-components 
network model for red cells ^1. If a curvature energy is introduced on the com- 
partment in place of the c ano nical bending energy S2 in the compartmentalized 
fluid surface model in Ref. we have also fluid surface models with cytoskele- 
tal structures ^3. A large variety of shape transformations are observed in such 
inhomogeneous fluid surface models, where the bond potential Si is the Gaussian 
bond potential, and the curvature energy S2 is the one-dimensional bending energy 
S^ 1 defined only on the compartments, which are one-dimensional objects linked 
with junctions . The phase structure depends on the elasticity at the junctions; 
a planar phase, and a tubular phase are observed in those models. The reason for 
such a variety of phases is closely connected to the cytoskeletal structure and the 
lateral diffusion of vertices. In fact, S^ 1 is considered to be insensitive to the surface 
shape, because S^ is defined only on the compartments in contrast to the model 
in Ref. where is defined all over the lattice. The surface shape is not always 
uniquely determined if the curvature is given only at small part of the surface, and 
moreo ver large surface fluctuations are expected in the compartmentalized model 
in Ref. ^ due to the lateral diffusion of vertices inside the compartments. 



2.2. Self- avoiding surface models 

The current studies that have been conducted on SA surfaces are considered to be 
still in the pioneer stage. In this subsection, we briefly comment on the existing SA 
surface models and the results of the numerical studies. The SA surface model is 
defined by a SA int eract ion, which is an extension of the Hamiltonian of the Edward 
model f or polymers The phase structu re of the SA models has been extensively 
studied l3jj3 2 | 33 | 34 | 35 | 36 | 37 | 38 | 39 | 4() | 4l El a i though the tota i num ber of studies are 

currently considered to be far smaller than those of the phantom surfaces. 

We have two types of SA models for numerical studies: the ball-spring (BS) 
model and the impenetrable plaquette (IP) model. The BS model is defined on two- 
dimensional networks, which are composed of vertices and bon ds con necting two 

101 IQOI 

nearest neighbor vertices by a hard-core and hard- wall potential 1 1 \ The size of 
ball as the vertices and the length of spring as bonds are constrained such that no 
vertex can move from one side of a triangle to the other side. This SA potential of 
the BS model is defined between all pairs of vertices, however, the simulations are 
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slightly less time-consuming than those of the IP model. The SA interaction of the 
IP model is defined such that the triangles are constrained to avoid intersecting. 
Although the simulations of the IP model are relatively time consuming, the IP 
model seems advantageous to the BS model. In fact, two neighboring triangles i 
and j of the IP model can completely bend such that 1— n; ■ n.; = 2, while in the BS 
model the bending angle 0y is constrained such that 9ij < 9q , where cos 6*^- = • n,, , 
and Qq{< tt) is determined by the SA potential. 

.■^fe. e _ cr umplirig; transition is reported to disappear from the SA FC surfaces 
\ 6 3 | 34 | 3o | 3o | 3 1 \ 3a | 44 is because no completely-collapsed phase appears in the 

SA FC s urfaces. In fact, no crumpled phase is observed in both of the BS model 
| 33 | 34 | 37 | and the ip model mMMM. 

To the contrary, the smooth phase is expected 
to remain unchanged from that of the phantom surfaces. However, numerical res ults 
are not always universal; in fact, the Hausdorff dimension of the IP model of ^Hl 
is H ~ 2.3, while that of Ref. is H = 2.1(1), although both models are defined 
on the disk surface. As mentioned in the Introduction, it is possible that the SA 
interaction plays a non-trivial role in membrane shapes in the smooth phase 
Thus, we should study the SA model more extensively. 

To summarize the comments including those in the previous subsection, we have 
several phantom surface models, which have a multitude of phases. The models 
are considered as non-trivial variations of the canonical surface model. The phase 
structures of almost all models have not yet been studied on the SA surfaces. The 
current understanding of the phase structure of SA surfaces are as follows: the 
crumpling transition disappears from the FC model, because the SA interaction 
prohibits the surfaces from collapsing in both of the BS model and the IP model. 
The smooth phase of the SA surfaces are considered to be almost identical to the 
smooth phase of the phantom surfaces, while the membrane shapes are expected to 
be influenced by the SA interaction under some specific conditions. 



3. Model 

In this section, we de fine a SA model, which corresponds to the phantom surface 
models in Refs. 129130] x r i an g U l a ted lattices of sphere topology are assumed to define 
the model, and the lattices are constructed using the icosahedron. By splitting the 
edges and faces of the icosahedron, we have a lattice of size N = l0£ 2 + 2, where £ 
is the devision number of an edge of the icosahedron. The coordination number q 
of vertices is q — 6 almost everywhere excluding 12 vertices of q = 5. The lattice is 
characterized by the three numbers N, N B (=3N-6 = 30£ 2 ), and N T {=2N^l = 20£ 2 ), 
which are the total number of vertices, the total numbe r o f bonds, and the total 
number of triangles, respectively. The lattices used in Ref. are random lattices, of 
which the coordination number is not always uniform and, they are slightly different 
from the lattices constructed as above. However, the phase structure of FC surface 
models is expected to be independent of the lattice structure ^S. 

The dynamical variable of the FC model is the position R 3 ) of the vertex 
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i(= 1, • • • ,N). The partition functions of the model are given by Eq. ([IJ. The 
Hamiltonian S(X) is defined by a linear combination of the area energy Si, a 
curvature energy S2, the pressure term — ApV, and a SA potential U, such that 

S(X) = Sx + aS 2 - ApV + U, 

S 1= J2 A ^ S 2 = -X>g(*/ 27r ). (3) 



U=Y,U(A,A'), U(A,A') 

A.A' 



00 (triangles A, A' intersect), 
(otherwise). 



Si is the sum over the area A a of triangle A. The symbol Si in S2 is the sum of 
internal angle of triangles meeting at the vertex i. S2 can be called a deficit angle 
energy, although S2 is different from the sum of the deficit angle 6{ — 2ir of the 
vertex i. If S2 is defined without "log" and is given by ^tt), then S 2 depends 

only on the surface topology and is a constant on piece-wise linearly triangulated 
surfaces. However, S 2 in Eq. (|3]) is well-defined as a curvature energy because of the 
log function as mentioned in Section [2Tl The symbol a[fcT] denotes the bending 
rigidity, where k is the Boltzmann constant and T is the temperature. 

V is the volume enclosed by the surface, and Ap is the pressure which is defined 
by Ap = Pi n — p ou t, where p ou t (Pin) is the pressure outside (inside) the surface. If 
p ou t is assumed to be p ou t = 0, then the positive (negative) Ap implies that p- ln is 
positive (negative). We should note also that the volume V is well defined only if 
the surface is self-avoiding. V is bounded below such that V > in the SA surfaces, 
while V can be negative in non SA surfaces. 




(a) "" (6) (c) 

Fig. 1. (a),(b) Two intersecting triangles, and (c) an intersection of the bond I and the triangles 
with the vertices X'-, Xj and X^, where X 1 - is a new position of the vertex i. In (a), the bonds 
AB and BC of the triangle ABC intersect with the triangle DEF, while no bond of the triangle 
DEF intersects with the triangle ABC. In (b), the bond BC of the triangle ABC intersects with 
the triangle DEF, and the bond DE of the triangle DEF intersects with the triangle ABC. 



A' m t ne SA potential U denotes the sum over all pairs of non nearest 
neighbor triangles A and A'. The potential U(A, A') is defined such that any pairs 
of non nearest neighbor triangles A and A' should not be intersecting. Figure [T^a) 
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shows two pairs of intersecting triangles, in which the triangle ABC penetrates 
the triangle DEF or in other words the bonds AB and BC intersect with the 
triangle DEF. On the contrary in Fig. [ljb), the triangle ABC and the triangle 
DEF intersects with each other, or in other words a bond of one triangle intersects 
with the other triangle and vise versa. We describe the numerical implementation 
of the SA interaction U in detail in the following section. 

The SA potential U in Eq. ^ is not identical to the one assumed in the SA model 
in Ref. ^21 gj-^ hence, the surface is completely self-avoiding under the potential 
U. In fact, the triangles are allowed to intersect with finite energy in Ref. while 
those in the model of this paper are prohibited to intersect with each other because 
of the infinite energy assumed in U. 

Finally in this section, we comment on how to compute the volume V enclosed 
by the surface. The initial value of V in the simulations is assumed such that 
V = 4-7rr 3 /3, where r is the radius of the initial configuration of sphere lattice. This 
initial value V = 4irr 3 /3 is slightly larger than the real volume, because the surface 
is linearly triangulated. However, it is almost evident that the deviation can be 
negligible in the limit of N — > oo. The volume V changes during the simulations 
according to the rule V — > V + AV every update of vertex, where AV is the volume 
of small tetrahedra, such as the one shown in Fig. QJc). AV is positive or negative, 
which is determined according to whether the new position X[ is outside or inside 
the surface, in which the orientation is uniquely fixed by a normal vector of each 
triangle. We should note that AV is well defined only when the surface is self- 
avoiding. It is apparent that AV is not well defined when some part of volume 
element of AV is shared by some other AV' , i.e., the surface is allowed to self 
intersect. 

The enclosed volume V can also be computed by using the divergence theorem 
applying the position vector of the center of mass of the triangle i. Not only 
AV but also V is exactly identical to the one obtained by the above mentioned 
technique. A very small deviation can be seen in the total volume V, however, it 
is less than 1% even in the cup like phase on the N = 1442 surface during the 
simulations. This small deviation of V is the one between y = 47rr 3 /3 and V of the 
initial triangulated sphere. 

4. Monte Carlo technique 

The canonical Metropolis Monte Carlo (MC) technique is employed for simulating 
the integrations of the variables X in of Eq. ([I]) . The three-dimensional random 
move X — > X' = X+SX is accepted with the probability Mm[l, exp(— SS)], where SS 
is given by SS = S'(new)— S'(old) under the constraint of the potential U. The symbol 
5X is randomly chosen in a small sphere, whose radius is fixed in the simulations 
such that the acceptance rate rx of X' should be approximately rx = 50%. 

The constraint of U(A, A') in Eq. Q is composed of two different constraints 
on a new vertex position as follows: let Xi and X[ denote the current position and 
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the new position of the vertex i as shown in Fig. [He). The shaded triangle in Fig. 
[TJc) forms a new surface. One of the constraint imposed on X[ is that the new 
triangle i'jk has no intersection with the disjoint bonds, where "disjoint bonds" are 
the edges of triangles disconnected with the triangle ijk. The other constraint is 
that every new bond, such as the bond i'j in Fig.[ljc), has no intersection with the 
disjoint triangles, where "disjoint triangles" are those disconnected with the bond 
ij. These two constraints imposed on X[ make the surface self-avoiding in the sense 
that any two disjoint triangles have no intersection with each other. 

The first constraint prohibits the new triangle i'jk shown in Fig. [Tic) from being 
penetrated by disjoint triangles. The second constraint imposed on X[ prohibits 
the new triangle i'jk from penetrating some other triangles. The intersection of 
the triangles shown in Fig. [Tfb) is prohibited by both of the constraints, while the 
intersection in Fig. HJa) is prohibited only by one constraint and is not prohibited 
by the other constraint. This is the reason why two constraints are necessary to 
make the surface self-avoiding by checking an intersection of a bond and a triangle. 

We assume a sphere of radius Rq at the center of mass of the triangle i'jk 
shown in Fig. [IJc), and check whether or not the triangle intersects with disjoint 
bonds inside the sphere. The check of intersection in the second constraint is also 
performed assuming the sphere of size Rq at the center of the bond i'j. The radius Rq 
is assumed to be Ro = 6(L) , where (L) is the mean bond length. As a consequence, 
the computational time is reduced by 20% ~ 60% or more, which depends on a. 

The bond length L and the triangle area are bounded below such that 
L > 1 x 10~ 7 and Aa > 0.5 x 10~ 7 in the simulations. The final results of the 
simulations are considered to be independent of these lower bounds, because these 
bounds are sufficiently small and almost all bond lengths and triangle areas are 
larger than these values. 

The total number of MC sweeps (MCS) after the thermalization MCS is about 
1 x 10 7 ~ 2x 10 7 on the N = 1442 surface, and relatively small number of MCS 
is assumed on the smaller surfaces. The total number of the thermalization MCS 
is about 0.5 x 10 6 . The thermalization MCS in the collapsed tubular phase is very 
large; it is sometimes 1 x 10 7 or more at the phase boundary close to the cup like 
phase on the N = 1442 surface. Intersection of bonds with triangles is checked 
every 5xl0 5 MCS throughout the simulation; the check is performed between every 
disjoint pair of bond and triangle. No intersection is observed at every assumed 
value of a including a = and Ap. 

5. Results 

The snapshots of FC surfaces at Ap — are shown in Figs. [2]Ja)-[2jf). The surface 
size is TV = 1442. The assumed bending rigidities are (a) a = 0, (b) a — 100, (c) 
a = 500, (d) a = 1000, (e) a = 1 x 10 4 , and (f) a = 2 x 10 4 . The scales of the figures 
are all different from each other. The surfaces shown in the figure are considered to 
be in (a),(b),(c) the wrinkled phase, (d),(e) the tubular phase, and (f) the smooth 



(a) (b) (c) (d) (e) ffi 




Fig. 2. The snapshots of FC surfaces and the surface sections of size N = 1442 obtained under 
Ap = at (a) a = (wrinkled), (b) a = 100 (wrinkled), (c) a = 500 (wrinkled), (d) a = 1000 
(tubular), (e) « = lxl0 4 (tubular), and (f) « = 2xl0 4 (smooth spherical). 

spherical phase. The surface in Fig. |2fa) can be called a collapsed surface because 
the surface is highly fluctuating, however, it encloses empty space inside the surface 
and, therefore, the surface is not always crumpled in the limit of a—}Q. The spherical 
surfaces in Figs. Hfb) and [2fc) look slightly smooth, however, they are apparently 
different from the surface at sufficiently large a shown in Fig. [UJf). The surfaces 
in Figs. [2Jd) and [^e) can be called a tubular surface. The surface in Fig. [^f) 
is very smooth and can be called the smooth spherical surface. All of the phases, 
excluding the wrinkled phase, correspond to those of the same model without the 
SA interaction in Ref. The collapsed phase can be seen in the model in Ref. I^S, 
while it is not in the SA model at least under Ap = Q. 

Snapshots of the FC surfaces and the surface sections are shown in Figs. [3]Ja)- 
|3ff), where a negative pressure Ap = —0.5 is assumed. The snapshots are slightly 
different from those at Ap = shown in Figs. [UJa)-[2jf). The snapshots in Figs.[3]Ja) 
and[3fb) indicate that the surfaces in the collapsed phase are almost crumpled. We 
see that the tubular surface in Fig. |3jc) is also collapsed. The cup like surfaces in 
Figs. E2d) andElJe) are new and typical of the condition Ap = —0.5, therefore, we 
call the new phase as the cup like phase. The smooth phase in Fig.|3jf) corresponds 
to the smooth phase in Fig. [2]Jf) at Ap = 0. The phase structure at a — > oo is 
understood to be independent of Ap. 

We also see that almost all parts of the surfaces in Figs. EJJd) and EJe) consist 
of oblong triangles and are locally smooth along one specific direction and wrinkled 
along the direction vertical to the smooth direction. This is also expected in the 
linear phase shown in Fig. EJc) at Ap = —0.5. To the contrary, the surface in the 
wrinkled phase shown in Fig. [2ja) consist of almost regular triangles and locally 
wrinkles along any directions. In the case of smooth phase in Fig. [2jf ) , the surface 
is smooth along any directions. Thus, the surface is symmetric under the three- 
dimensional rotations both in the limit of a —¥ oo and a — > 0, while the rotational 
symmetry is spontaneously broken at intermediate region of a. This observation 
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Fig. 3. The snapshots of FC surfaces and the surface sections of size N = 1442 obtained under 
Zip = —0.5 at (a) a = (collapsed), (b) a = 100 (collapsed), (c) a = 1.5xl0 4 (collapsed tubular), 
(d) a = 2xl0 4 (cup like), (e) o = 3xl0 5 (cup like), and (f) a = 4 X 10 s (smooth spherical). 



is independent of the two values of Ap. This symmetry breakdown or restoration 
is considered to be closely connected to the structural change of the constituent 
triangles; the symmetric surfaces are composed of almost regular triangles, while 
the non-symmetric surfaces are composed of oblong triangles. 
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Fig. 4. [Si - (3/2) Ap V]/N vs. a under (a) Ap = and (b) Ap = -0.5. The error bars on the 
symbols denote the standard errors. The solid lines connecting the symbols are drawn as a guide 
to the eyes. 
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In the following presentations, we show how the shape transformation transitions 
and/or the SA interaction are reflected in the physical quantities including the 
Hausdorff dimension H in the limit of a— >0. 

First of all, we show [Si-{3/2)ApV]/N, denoted by S[/N, in Figs. Ha) and 
IDJb). Because of the scale invariant property of the partition function Z& x of Eq. 
(HJ), S[/N is expected to be S[/N = 3/2 at sufficiently large TV. We see that all of the 
results are consistent with the prediction. This implies that the volume V is well- 
defined and that the SA interaction is correctly implemented in the simulations. We 
note that it is straightforward to prove that S[/N — 3/2 4 As described in Section 
12.11 the scale invariance of Z is represented by dZ(XX)/d\\\=i = 0. Because of the 
scale transformation X — > XX, Si and V change to A 2 Si and A 3 V^ while S 2 and U 
remain unchanged. Since the integration J ]X dXi also changes to A 3 ^ -1 ) J Y[i dXi, 
then we have the relation S[/N = 3/2 in the limit of TV— s-oo. 

S[ is identical with 5*1, which is the total area of surface, in the case Ap = 0, and 
therefore, the scale invariance implies that the surface area remains unchanged in 
the whole range of a. To the contrary, the surface area Si discontinuously changes 
at the transition points in the case Ap = — 0.5 at least, because V discontinuously 
changes at the transitions as we will see below, while S[ remains unchanged. This 
implies that the internal property of surface is significantly influenced by the exter- 
nal condition Ap. 
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Fig. 5. The volume V vs. a under (a) Ap = and (b) Ap = — 0.5. The vertical dashed lines denote 
the phase boundaries of the N = 1442 surface. The solid lines connecting the symbols are drawn 
as a guide to the eyes. 



The volume V enclosed by the surface should be bounded below such that V > 0, 
which is satisfied only if the surface is self-avoiding. The model in this paper is 
strictly self-avoiding, and hence V is expected to be well defined even when Ap is 
large negative. Figures [SJa) and[SJb) show the dependence of V on a under Ap = 
and Zip = —0.5, respectively. The vertical dashed lines in the figures represent the 
phase boundaries between two different phases just like in Fig. |4] The name of the 
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phases corresponds to the surface shape, which can be visualized as snapshots just 
like those in Figs. [2] and [3] 

The detailed informations such as the order of the transitions are not obtained. 
It is possible to perform the finite-size scaling analyses to see the order of the 
transitions by performing the simulations at the transition region more extensively 
however, we confine ourselves of the phase structure in the wide range of a and, as 
a consequence, the order of the transitions is not fully examined. Thus, it remains 
unclear whether or not the smooth spherical phase and the tubular phase (or the 
cup like phase) are separated by a first-order transition, although the volume V 
discontinuously changes at the phase boundary. We see that the volume V in the 
collapsed phase is larger than that in the cup like phase under Ap =—0.5 at least, 
while V at a— s-0 is smaller than that in the tubular phase under Ap = 0. 




Fig. 6. The mean square size X 2 vs. a. The vertical dashed lines denote the phase boundaries of 
the iV = 1442 surface. 



The mean square size X 2 is defined by 

i i 

where X is the center of mass of the surface. The value of X 2 changes depending 
on the distribution of the vertices in R 3 , and hence X 2 as well as V can reflect 
shape transformations. However, the quantity X 2 does not always show the same 
behavior against a as that of V. Figures [6ja) andEJb) show X 2 vs. a under Ap = 
and Ap = —0.5. We see in Fig. [6ja) that X 2 discontinuously changes at the phase 
boundary between the smooth spherical phase and the tubular phase. It is also 
easy to see from Fig. E|b) that X 2 discontinuously changes at the phase boundaries 
between the smooth spherical phase, the cup like phase, and the collapsed tubular 
phase. 

Figures EJa) andJT^b) show X 2 obtained under Ap = Q and Zip = —0.5 at small 
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Fig. 7. The mean square size X 2 vs. a at small a region under (a) Ap = and (b) Ap = —0.5, 
and (c) X 2 vs. TV in a log-log scale obtained at a = under Ap = and Zip = —0.5. The straight 
lines in (c) are drawn by fitting the largest three data points to Eq.lBl. 



a = a = 0.1 a = 1 a = 2 

Zip = if = 2.33±0.08 if = 2.34±0.08 if = 2.33±0.08 
Z\p=-0.5 if = 2. 60±0.17 - - if = 2. 59±0. 17 



a region. The Hausdorff dimension if of the surface is defined by 

X 2 ~N 2 / H (N-too). (5) 

By fitting the data X 2 obtained at a = to Eq. we draw straight lines in Fig. 
|7fc), and the values of if are shown in Table Q] including if at a = 0.1 and a = 1 
under Zip = 0, and at a = 2 under Ap = —0.5. The fitting is performed using the 
largest three data points under each condition of Zip. We have if ~ 2.6 at Ap= —0.5 
and H ~ 2.33 at Z\p = in the limit of a — > 0. The value of if ~ 2.33 is compatible 
with the one if ~ 2.3 in Ref. l^SI, while if ~ 2.6 at Zip = —0.5 is slightly larger 
than the Flory estimate 2.5 and compatible with the fact that the surface is almost 
crumpled as we see in the snapshot in Fig. Efa). We should note that the value 
H ~ 2.3 seems independent of the details of the model, the SA interaction, and the 
surface topology. However, the result if ~ 2.33 is larger than the one if = 2.1(1) of 
Ref. SH, thus it is also possible that H depends on the model on the SA surfaces. 

We comment on the size effect of the results in Table Q] As mentioned above, 
the data obtained on the small sized surfaces, such as N = 162 and iV = 362, were 
excluded from the fitting. By including the small two data in the fitting, we have 
H = 2.50(5) for Zip = 0, a = and if = 2.93(7) for Zip = -0.5, a = 0. Both of if are 
slightly larger than H = 2.33(8) and ff = 2.60(17) shown in Table [TJ Thus, the size 
effect is not negligible at least on the surfaces iV<362. 

The surface size can also be reflected in the maximum linear extension Le, which 
is defined by the maximum distance between two vertices on the surface: 



L E = Mnx{\X t -X J \ | = !,••• ,N)}, 



(6) 
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Fig. 8. The maximum linear extension Le vs. a. The dashed lines denote the phase boundaries 
of the ^ = 1442 surface. 



where Xi and Xj are not always connected by a bond. The phase transition of shape 
transformation is also reflected in the structure of triangles; we see in the snapshots 
in Figs. [2] and [3] that the surface consists of almost regular triangles in the smooth 
spherical phase while it includes oblong triangles in the tubular phase, where Le 
expected to be very large. We expect that this structural change is reflected in Le- 
FigureslHJa) and[8jb) show Le vs. a under Ap = and Ap=— 0.5. The discontinuous 
change of Le at the phase boundaries shown in the figure implies that the phase 
transitions are accompanied by the structural cha nge of sur faces . This structural 
change is typical of the Nambu-Goto surface model B USSED] 
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Fig. 9. The two-dimensional bending energy Ss/Nb vs. a under (a) Ap = and (b) Ap = — 0.5. 



The two-dimensional bending energy S3/NB is shown in Figs. IDJa) and[9]^b), 
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where S3 is defined by using a unit normal vector of the triangle i such that 



We write the two-dimensional bending energy as S3 to distinguish it with the deficit 
angle energy S2 in Eq. ([3]). 

We see that S^/Nb discontinuously changes at the phase boundaries, where the 
physical quantities such as V, X 2 and Le discontinuously change. To the contrary, 
the deficit angle energy S2 defined in Eq. ([3]), which is not shown in the figures, 
appears to vary almost smoothly in the whole range of a. At the boundary between 
the smooth spherical phase and its neighboring phase, S2/N is expected to change 
discontinuously like the other physical quantities. However, the discontinuity is very 
small and it is almost invisible just as in the case of the self-intersecting model in 
RefM. 

6. Summary and Conclusion 

We have numerically studied a self-avoiding (SA) surface model on fixed- 
connectivity (FC) triangulated lattices of sphere topology. The self-avoidance of 
the model in this paper is not identical to those of the well-known SA models; the 
ball spring model and the impenetrable plaquette (IP) model. However, the SA 
model in this paper belongs to the IP models, because the intersection of disjoint 
triangles are prohibited by the SA interaction. The phase structure of the FC model 
under Ap = is found to be almost identical to that of the phantom surface model 
in Ref. 1221 except for the evidence that the collapsed phase disappears from the SA 
model. Thus, the influence of the SA interaction on the phase structure is very small 
contrary to the expectation that the SA interaction can suppresses the multitude 
of phase transitions in the phantom surface model. 

To be more precise, the model in this paper is a Nambu-Goto surface model with 
a deficit angle energy. The SA interaction is defined such that all possible pairs of 
non-nearest neighbor triangles are prohibited from intersecting. Because the volume 
enclosed by the SA surface is well defined, the pressure term — Ap V can be included 
in the Hamiltonian. The simulations are performed under Ap = 0, and Ap =—0.5 
on the FC surfaces, where Ap=— 0.5 implies that the pressure inside the surface is 
lower than the pressure outside the surface. 

Our observations on the FC surfaces are as follows: the smooth spherical phase, 
the tubular phase, and the collapsed phase can be seen under those two conditions 
of Ap, and the cup like phase is seen under zip = —0.5. Thus, the phase structure of 
the model under Ap = is almost identical to that of the phantom surface model, 
although the collapsed phase is slightly different from each other; the collapsed 
surfaces are completely shrunk in the phantom surface model, while the SA surfaces 
are not completely shrunk at Ap = at least. The Hausdorff dimension H = 2.33(8), 
obtained at a — under Ap — 0, is independent of the curvature energy and is 
considered as the Hausdorff dimension of the Nambu-Goto SA surface. This result 




(7) 
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is consistent with the known result of H ~ 2.3 of the IP model in Ref. ^3H1, where 
the model, the SA interaction, and the surface topology are different from those in 
this paper. In this sense, it is possible that the value H~2.3 depe nds only on the 

IOQI 

self- avoidance, although the surface size of the simulation in Ref. 1:221 is relatively 
smaller than those assumed in this paper. To the contrary, H = 2.33(8) is larger 
than the result H = 2.1(1) of Ref. S3, and therefore, it is also possible that H of the 
SA surface depends on the model. The SA surface models should be studied more 
extensively. 

It is also interesting to study whether or not the mul titude of phases in the 
fluid surface models with cytoskeletal structures in Ref. ^ is observed under a SA 
interaction. The SA interaction assumed in the model of this paper can also be 
assumed in those fluid surface models even when Ap is negative. This remains to 
be a future study. 
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